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ABSTRACT 

Day and night mapping of the global distributions of 
the horizontal cloud-cover and the corresponding cloud- 
top pressure levels can be derived from the same infra- 
red data used to derive clear column temperature profiles. 
Application to the 15 urn VTPR data are given. Extension 
of this approach for the determination of the radiative 
transfer properties of clouds are presented and the pos- 
sibility of using such information to infer cloud types 
is discussed. 

1 . INTRODUCTION 

The most important atmospheric parameters of interest in 
studying radiative transfer in cloudy atmospheres are the dear- 
column vertical temperature profiles, the amounts, heights, and 
radiative transfer properties of clouds. Clouds play a major role 
in the radiative processes in the earth atmospheres both in the 
absorption and reflection of solar radiation and in the emission 
of thermal energy. A knowledge of the vertical location, the hori- 
zontal distribution and the optical properties of clouds provides 
information on the heat sources and sinks, storage rates and trans- 
port phenomena in the atmosphere. Such%information play a critical 
role in determining the drives for motions in the atmosphere and 
oceans. 

Visible and infrared pictures of clouds indicate the exis- 
tence of clouds, but do not give information about their altitudes 
or their radiative transfer parameters which are essential for the 
estimate of their types. By contrast, infrared temperature sound- 
ing data permit the determination of the horizontal and vertical 
distribution of clouds provided that all the sounding channels ob- 
serve the same field of view at the same time. The requirement of 
simultaneous observations is necessary because clouds are usually 
a very inhomogeneous medium both vertically and horizontally. 

2. APPROACH 

The upwelling radiance from a planetary atmosphere is a 
function of the thermal state of the atmosphere, the concentration 
of radiatively active gases, and the extents, heights, and radiative 


197 


transfer properties of clouds and aerosols. Thus, in principle, 
it should be possible to recover useful information about the 
physical and chemical structure of an atmosphere from analysis of 
the upwelling radiance. However, the problem in analyzing such 
data lies in finding ways to uncouple the effects of these vari- 
ables and retrieve the true values of each unknown separately, oy 
U&cuUng thz cloud cUccZi cu> ikofit tcnm o&ciMationi oveA the 
cleoA-column upoMing mdiance, an analytical method was devel- 
oped by Chahine (1974) to retrieve dear-column vertical tempera- 
ture profiles from radiance measurements made in the presence of 
clouds. The method requires radiance data measured in two spec- 
tral regions and over two adjacent fields of view having different 
amounts of clouds. The uncoupling of the effects of clouds is 
carried out without any a priori information about the amounts, 
heights and optical properties of the clouds in the fields of view. 
Once the dear-column temperature profiles are determined the same 
radiance data could then be used to determine the heights, amounts, 
and radiative transfer properties of clouds. 

Formulation of this approach is straightforward. It will be 
carried out here in two steps: 

First, we consider the radiance I(v) measured at frequency 
V in the presence of clouds, and define the corresponding dear- 
column radiance I(v) as 


I(v) = B(v,T5) t(v,Pj) + r B[v,T(p) In p ( 1 ) 

in P3 

where t(v,p) is the transmittance of a dear-column of gaseous ab- 
sorbers between pressure level p and the sounder. B is the black- 
body Planck function and T(p) is the dear-column vertical tempera- 
ture profiles. Next we express the difference between I(v) and 
I(v) in terms of an expansion function G(v,p) and an expansion co- 
efficient N as 

I(v) - I(v) = N 6 (v,p^, ...) + N'G'(v,p^ ...) (2) 

In this first step we don't need to define the form of G(v,p) be- 
cause we aim to eliminate it. Detailed discussion of this approach 
is given by Chahine (1977). For simplicity we take here one ex- 
pansion team only and consider observations made over two adjacent 
fields of view (subscripts 1 and 2) having different amounts of 
clouds, such that Ii f I2 and write Eq. (2) as 

Ii('^) - Ii(v) = Nj G(v,p, ...) 

l2(v) - l2(^^) “ ^2 G(v,p, ...) 

Now, if the two fields of view are continguous we can assume that 
Ii(v) ^ l2(^^) ^nd drop their subscripts. Equation (3) becomes then 
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( 4 ) 


I(v) = Ii(v) + n[Ii(v) - l2(v)] 

where n = „ „ = constant 

N2 - Nj 

Determination of T(p) and I(v) could be carried out according to 
the methods described by Chahine (1974, 1975, 1977). 

In the second step we substitute the value of I(v) into 
Eq. (3) and write 


Ni 6(v,p^, ...) = I(v) - Ii(v) (5) 

for each sounding frequency vi. The left-hand side of Eq. (5) con- 
tains all the radiative transfer properties of the clouds and the 
terms on the right-hand side are all known. We need now to formu- 
late or model the expansion function G(v,p ...). 

3. CLOUD MODELING 

To deduce the properties of clouds one is faced with the 
problem of solving the complete equation of transfer in a Mie scat- 
tering medium. This sort of approach is not practical here because 
of the small number of available measurements. Thus we are con- 
fined to use simple black or nonblack cloud models consisting of 
one, two or three layers. 

For the types of data currently available from infrared 
sounders observations made in the range between 3.7 pm and 15 pm, 
we need consider only one layer cloud models. We could express 
G(v,p, ...) in terms of only a small set of cloud parameters such 
as the fractional cloud cover N^., the mean cloud top pressure level 
P(-, the cloud emissivity e^., transmissivity t^., and reflectivity 
Pc- 


a. Single Layer of Black Clouds 


In most radiation models it is assumed for simplicity pur- 
poses that all clouds, except cirrus, are sufficiently thick and 
dense to be treated as black clouds with emissivity 6^=1. Low 
stratus and cumulus water clouds have emissivities which are close 
to unity. Yamamoto et al . (1970) have shown that for long-wave 
radiations dense water cloud absorbs more than 90% of infrared 
radiation within a depth of 50 m. Thus if we set the source func- 
tion at the cloud level P(~ to be the Planck function value B[v, 
Tc(Pc)] we can express G(v,p) for e^. = 1 as 


G(v,p^) = B(v,T^) t(v,p^) 


J in p 

^B[v,T(p) ^^i^d In p (6) 

in p. 


- B[v,T(p^)] t(v,p^) , 
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where the only unknown on the right-hand side is the cloud top 
pressure level Pr. 

Determination of p^ can be accomplished by considering two 
sounding frequencies, say and V 2 for which x(v,Ps) is prefer- 
ably small in order to minimize the effect of changes in earth sur- 
face emissivity. From Eq. (5) and Eq. (6) we write, for, say, the 
first field of view 


G(vi,p^) I(vi) - Ii(vi) 

I(v 2) - Ii(v2) 


and solve for pr by trial and error in a manner similar to the ap- 
proach described in Section 4.b by Chahine (1975). 

The values of the expansion coefficient N in this case are 
the effective fractional cloud cover which is the product of the 
actual fractional cloud cover and the actual cloud emissivity, 

N N can be readily derived from Eq. (5) using vj or ^2 

as 



Application of this model to the analysis of the VTPR sounder data 
has been carried out and the results will be discussed in Section 4. 


b. Single Layer of Transmitting Clouds 


The assumption of black clouds is attractively simple but 
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of stratocumulus clouds do not become black until their thickness 
exceeds half a kilometer. The emissivity of middle-level and cir- 
rus clouds remains less than one even when their thickness exceeds 
several kilometers. A useful procedure to treat non-black clouds 
is to assume the cloud to be transmitting but not reflecting with 


e^(v) + t^(v) = 1 

By considering observations in two spectral regions, say [v] = 4.3 
urn and [v']= 15 pm CO 2 bands, we can derive an expression for the 
ratios of emissivities in the two bands as 

e^,(v) I(v) - Ii(v) G(v',p^) 

I(v') -I.(v') 

Since the cloud emissivity is frequently dependent and a strong 
function of the water content of clouds and their thickness, the 
ratios of the emissivity can be used for deducing the cloud opti- 
cal thickness and identify their types. The determination of N, 
e(v)/e(v') and p^ is carried out simultaneously by the method 
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of trial and error starting with p^, <. Pg in a manner slightly more 
elaborate than the one followed for black clouds. Numerical simu- 
lations of this approach are currently under way and the results 
will be applied to the HIRS II data. 

4. APPLICATION OF BLACK CLOUD MODEL TO VTPR DATA 

Satellite verification of this approach has been carried out 
and the results are shown in Figs.- 1, 2, and 3. The computations 
were carried out, in collaboration with M. Halem, J. Susskind, and 
J. Forkash, using 15 pm VTPR data from the NOAA 4 satellite. The 
determination of the clear column radiance for the VTPR (i.e. the 
clear column temperature profile) requires the use of a priori as- 
sumptions about the surface temperature. The surface temperature 
here was obtained from the NOAA Analysis. The effects of errors 
in Ts on the retrieved values of pc is small particularly for pc 
< 700 m according to McCleese (1976). 

The resulting sample of cloud distribution shown here corre- 
sponds to one week of VTPR data from January 1 - 7, 1975 averaged 
on the GSFC-GLAS model grid size of 4° latitude by 5° longitude. 

The effective fractional cloud cover distribution varied 
between 0% and 89%. The graphical illustrations in Fig. 1 show 
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Fig. 1 — Map of effective infrared fractional cloud cover 
derived from the VTPR Sounder data for the period 
of January 1 - 7, 1975. 
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this distribution in four shades corresponding to the ranges less 
than 20%, 20% -50%, 50% -80%, and greater than 80%. The results 
obtained by this technique agree very well with the results ob- 
tained by J. C. Sadler of the University of Hawaii for the same 
period of time. 

The cloud top pressure levels varied between surface and 
145 mb and the results illustrated in Fig. 2 show the cloud top 


I I )000 MB - 500 MB |H 500 MB - 230 MB HH > 250 MB 
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Fig. 2 — Map of effective cloud top pressure level dis- 
tribution derived from VTPR Sounder Data for the 
period of January 1 - 7, 1975. 

level distribution in three shades corresponding to surf ace- 500 mb, 
500 mb- 250 mb and above 250 mb. These results possess the basic 
features of the global circulation and compare well with the small 
amount of available information. Additional verifications are 
still required. 

Statistical averages are presented in Fig. 3 which give the 
meridional profiles of zonally averaged cloud distribution, in 
fractions, and the cloud top pressure levels in mb: The Global 
average of the effective cloud cover as seen in the 15 ym band is 
shown to be 0.41, corresponding to 0.40 for the northern hemi- 
sphere and 0.42 for the southern hemisphere. 

The standard deviation of the fractional cloud cover for 
the 4° wide zone is given also in Fig. 3. The average zonal stan- 
dard deviation is 0.24 which is close to the value corresponding 
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Fig. 3 — Meridional profiles of zonally averaged infrared 
cloudiness for January 1-7,J975. 

to a constant probability distribution p(X) = 1 on [0,1]. By con- 
trast, the standard deviation of the meridionally averaged frac- 
tional cloud distribution is 0.1. Thus, the dispersion (variance) 
of the meridional distribution of clouds is six times larger than 
the zonal dispersion. 

5. CONCLUSION 

The results shown in Figs. 1, 2, and 3 demonstrate the capa- 
bility of temperature sounding radiance data to yield meaningful 
information about the cloud distribution even under black cloud 
model assumptions. More extensive verifications still need to be 
carried out. 

Extension of this technique to transmitting clouds will pro- 
vide information about the true cloud distribution and emissivity. 
The emissivity could be used in conjunction with the cloud height 
and amount of yield information about the types of clouds. 

Further extension of this technique to the cases of reflect- 
ing clouds and to multiple cloud layers is, in principle, possible. 
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but application of these cloud models to current sounder data is 
not meaningful because of the poor vertical resolution provided 
by current operational sounders. 
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